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ABSTRACT 

We report the discovery of a Jupiter-size planet transiting a relatively bright {V = 11.56) and 
metal-rich early K dwarf star with a period of ~ 2.9 days. On the basis of follow-up photometry 
and spectroscopy we determine the mass and radius of the planet, HAT-P-3b, to be Mp = 0.599 ± 
0.026 Mjup and Rp = 0.890 ± 0.046 Rjup. The relatively small size of the object for its mass implies 
the presence of about 75 worth of heavy elements (~ ^ of the total mass) based on current theories 
of irradiated extrasolar giant planets, similar to the mass of the core inferred for the transiting planet 
HD 149026b. The bulk density of HAT-P-3b is found to be pp = 1.06 ± 0.17 g cm'^, and the 
planet orbits the star at a distance of 0.03894 AU. Ephemerides for the transit centers are Tc — 
2,454,218.7594 ± 0.0029 -f iV X (2.899703 ± 0.000054) (HJD). 

Subject headings: binaries: spectroscopic — planetary systems — stars: fundamental parameters — 
stars: individual (HAT-P-3) — techniques: spectroscopic 

1. INTRODUCTION 

Ground-based surveys for extrasolar transiting planets 
have produced 18 discoveries to date. This sample is still 
small enough that individual discoveries often advance 
our understanding of these objects significantly by push- 
ing the limits of parameter space, either in planet mass, 
radius, or some other property. Here we announce a new 
transiting planet, the third to come out of the HATNet 
survey, around a star previously known as CSC 03466- 
00819. It is the smallest yet discovered photometrically, 
and appears to have a heavy-element content represent- 
ing a substantial fraction of its mass. 

2. PHOTOMETRIC DETECTION 

GSC 03466-00819 is located in our internally labeled 
field "G145" centered at a = 13'^48'", S = -h45°00', 
which was observed using the 11-cm aper ture HAT-5 
instrument of the HAT network (HATNet; iBakos et all 
120021 12004f ). located at the F. L. Whipple Observatory 
(FLWO) on Mt. Hopkins (Arizona). The field was ob- 
served in single-station mode, i.e., without networked op- 
erations between Arizona and HATNet 's other station 
in Hawaii. Observations over an extended interval of 
six months (from 2006 January 28 through 2006 July 
22) yielded some 3200 photometric measurements per 
star for approximately 10,000 stars in the field down to 
/ ^ 14. Following standard frame calibration procedures 
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Fig. 1. — (a) Unbinned instrumental /-band discovery light 
curve of HAT-P-3 obtained with HATNet, folded with the period 
of P = 2.899703 days, (b) Unbinned instrumental Sloan i-band 
photometry collected on UT 2007 April 27 with the 1.2-m tele- 
scope at FLWO, along with our best- fit transit model curve (see 
text). 



astrom etry was performed as described in iPal &: BakosI 
()2006l ). and aperture photometry results were subjected 
to the de-trending algorithms we refer to as Exter- 
nal Parameter De correlation (EPD, described briefly in 
Bakos et al.1l2007D. and Trend Filtering Algorithm (TEA; 
Kovacs et al.l I2OO5I) . We searched the light curves for 



box-shaped t r ansit signals using the BLS algorithm of 
iKovacs et al.l (|200l). A prominent 2.9-day signal was 
detected in the / « 10.52 magnitude star GSC 03466- 
00819 (also known as 2MASS 13442258-^4801432; a = 
13'^44™22';58, S = -h48°01'43'.'2; J2000) with a tran- 
sit depth of 0.01 4 mag a nd a dip- signif icance parameter 
(DSP) of - 11 (|Kovacs fc Bakof^H^. This discovery 
light curve is shown in Figure (TJi. The standard devia- 
tion of the measurements after the EPD procedure was 
0.009 mag. Altogether 5 partial and 2 full transits were 
observed. As is shown in the following sections the tran- 
siting companion is a Jovian planet. Hereafter we refer 
to the star as HAT-P-3, and to the planetary compan- 
ion as HAT-P-3b. The near infrared color of the star 
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{J — Kg ~ 0.49) corresponds to an early K dwarf, and 
the annual proper motion listed in the Tycho-2 catalog 
(fic^cosS = -23.3 ± 2.9 m as yr-\ ^ = -24.0 ± 2.6 
mas vr~^: |H0g et al.ll2000[ ) is consistent with this. The 
photometric ephemeris was determined by combining all 
recorded HATNet transit events with a full transit ob- 
served more recently at higher precision with a larger 
telescope, as described below. The period we obtain is 
P = 2.899703 ± 0.000054 days, and the reference epoch 
of transit center is = 2,454,218.7594 ± 0.0029 (HJD). 
We use this ephemeris for the remainder of the paper. 

We note that HAT-P-3 has a faint (/ w 15) and red 
(J — Kg — 0.86) companion at an angular separation 
of -10" (2MASS 13442345+4801387), which is entirely 
within the aperture used for the HATNet photometry. 
If this happened to be an eclipsing binary, we estimate, 
based on the relative fluxes, that 1.6-mag deep echpses 
of this star could mimic the observed shallow transits 
seen in the merged light curve. However, our follow- 
up photometry described in Sj3]with much better spatial 
resolution shows that the shallow dips in brightness are 
indeed present in the light curve of HAT-P-3, whereas 
the faint companion is constant at the 0.05 mag level. 

3. FOLLOW-UP OBSERVATIONS 

HAT-P-3 was observed spectroscopically with the 
FLWO 1.5-m Tillinghast reflector in order to rule out 
the possibility that the observed drop in brightness 
is caused by a stellar companion in an eclipsing bi- 
nary. Three spec tra were obtai ned with the CfA Dig- 
ital Speedometer (jLathamlll992l ) over an interval of 25 
days. These observations cover 45 A in a single echelle 
order centered at 5187 A, and have a resolving power of 
A/AA « 35,000. Radial velocities were obtained by cross- 
correlation and have a typical precision of 0.4 kms~^. 
They showed no variation within the uncertainties, rul- 
ing out a companion of stellar mass. The mean velocity 
is -23.8 ±0.1 kms"\ 

Higher resolution spectroscopy was conducted with the 
HIRES instrument (|Vogt et al.l(l994D on the Keck I tele- 
scope. With a spectrometer slit of 0'.'86 the resolving 
power is A/AA « 55,000, and the wavelength coverage is 
~3800-8000 A. An iodine gas absorption cell was used 
to superimpose a dense forest of I2 lines on the stel- 
lar spectrum and establish a highly accurate wavelength 
fiducial (see Marcy & Butler 1992). A total of nine 12- 
minute exposures were obtained between 2007 March and 
June with the iodine cell, along with one without I2 for 
use as a template. Relative radial velocities in the Solar 
System ba rycentric frame were derived as described by 
iButler et~a l. (1996), including full modeling of the spatial 
and temporal variations of the instrumental profile. The 
results, listed in Table [U have internal errors < 3 ms~^. 

The iodine-free template spectrum from Keck was used 
to infer the atmospheric parameters of the star. A 
spectral synthe sis modeling was ca,rried out using the 
SME software (jValenti fc Piskunovl fT99i l. with wave- 
length ranges an d atomic line data as described by 
IValenti fc Fischeil (|2005f ). The effective temperature 
(Teff), surface gravity (logg*), iron abundance ([Fe/H]), 
and projected rotational velocity (wsini) we obtain are 
listed in Table [21 where the uncertainties represent inter- 
nal errors. The temperature and surface gravity corre- 



TABLE 1 

Relative radial velocity measurements of HAT-P-3. 



BJD RV o-Rv 

(2,400,000+) (ms-i) (ms"^) Phase 



54187.01101 -38.9 1.5 0.051 

54187.12931 -67.4 1.6 0.092 

54187.97909 -78.9 1.8 0.385 

54188.05795 -59.3 1.3 0.412 

54188.98084 -1-75.9 1.4 0.730 

54189.06253 +75.9 1.4 0.759 

54189.11203 +73.9 1.5 0.776 

54250.85697 -48.1 2.1 0.069 

54278.86111 +65.2 2.9 0.727 



TABLE 2 
Stellar parameters for HAT-P-3. 



Parameter Value Source 



Teff (K) 5185 ± 46 SME'' 

[Fe/H] +0.27 ± 0.04 SME 

logg* (cgs) 4.61 ± 0.05 SME 

iisini (kms-l) . . . . 0.5 + 0.5 SME 

(M0) 0.936tg:g^^ Y^+LC+SMEb 

R^, (Kq) 0.824^°;°^^ Y2+LC+SME 

(Lq) 0.442t°;g^^ Y^+LC+SME 

Mv (mag) 5.86 ± 0.20 Y^+LC+SME 

Age (Gyr) 0.4lgj Y^+LC+SME 

Distance (pc) 140 + 13 Y^-t-LC+SME 



^ SME = 'Spectroscopy Made Easy' packa ge for analysis 
of high-resolution spectra IValenti &: Piskunovl l|1996l ). See 
text.!^ Y^+LC+SME = Yale-Yonsei isochrones l IYi et al.ll200ll) . 
light curve parameters, and SME results. 

spond to an early K dwarf. 
Photometric follow-up of HAT-P-3 was carried out in 

the Sloan i band with KeplerCam (see Holman et al.l 
l2007t ) on the 1.2-m telescope at FLWO. After several 
attempts thwarted by weather we obtained a reasonably 
complete and high-precision light curve of HAT-P-3 on 
UT 2007 April 27. An astrometric solution between the 
individual frames and the 2MASS catalog was carried 
out using first order polynomials based on ~100 stars 
per frame. Aperture photometry was performed using a 
series of apertures on fixed positions around the 2MASS- 
based [AT,!"] pixel coordinates. We selected a frame taken 
near the meridian and used ~70 stars and their magni- 
tudes as measured on this reference frame to transform 
all other frames to a common instrumental magnitude 
system (weighting by the estimated photometric errors 
of each star). The aperture yielding the lowest scatter 
outside of transit was used in the subsequent analysis. 
The light curve of HAT-P-3 (Figure [T^) was then de- 
correlated against trends using the out-of-transit sections 
and a dependence on hour angle and zenith distance. 
The photometry is somewhat affected by occasional thin 
cirrus during the night, and has a typical precision of 
~2.5 mmag per individual measurement at 20-second ca- 
dence. 

4. ANALYSIS 

Modeling of the z- band light curve was ca rried out us- 
ing the formalism of 'M andel fc Agol (2002). Quadratic 
limb-darkening coefficients (ui = 0. 3832 and — 
0.2683) were taken from the tables of lClaretl ()2004f) by 
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TABLE 3 

Spectroscopic and light curve solutions for HAT-P-3, and 
inferred planet parameters. 

Parameter Value 

Spectroscopic parameters 

P (days)'' 2.899703 ± 0.000054 

Tc (^0-2,400,000)=^ 54,218.7594 ± 0.0029 

K (ms-l) 89.1 ± 2.0 

7 (ms-i) -14.8 ± 1.5 

e (adopted) 

Light curve parameters 

a/R. 10.5910-^^ 

Rp/R* o.iioolHSS 

b = acosi/Ri, O.Slto^g 

i (deg) 87.24 ± 6.69 

Transit duration (days) 0.0858 ± 0.0020 

Planet parameters 

Mp (Mjup) 0.599 ± 0.026 

Rp (R.iup) 0.890 ± 0.046 

Pp (g cm-=*) 1.06 ± 0.17 

a (AU) 0.03894 ± 0.00070 

loggp (cgs) 3.310 ± 0.066 

Held fixed from the photometric determination (i[2t. 

interpolation to the stellar properties given in Table O 
The period was held fixed, and the orbit was assumed to 
be circular, since the eccentricity is not sufficiently well 
constrained by the presently available velocity measure- 
ments. The four adjusted parameters are the planet-to- 
star ratio of the radii (i?p/i?*), the normalized separation 
{a/ Ri,) where a is the semimajor axis of the relative or- 
bit, the impact parameter {h = acosj/i?*), and the time 
of the center of the transit (Tc). A grid search was per- 
formed at a fixed to locate the minimum in the 
parameter space of the other three variables, and these 
were subsequently held fixed to optimize T^.. The process 
was iterated until convergence. Uncertainties were esti- 
mated by alternately fixing each parameter at a range 
of values near the minimum and allowing the others to 
float, until the best- fit solution produced an increase in 
the corresponding to a Icr change. The results are 
shown in Tabled 

The relative radial velocities of HAT-P-3 were fitted 
with a Keplerian model constrained to have zero eccen- 
tricity. The ephemeris was held fixed at the value for 
P determined in Sj^land from the light curve analysis. 
The solution fits the data well (see Figure[2^), but has an 
rms scatter of 4.7 m s~^ (accounting for the number of de- 
grees of freedom) that is larger than the typical internal 
errors (1.3-2.9 ms~^). No obvious long-term trends are 
seen in the residuals. Inspection of the spectra reveals 
that the Ca II H and K lines display modest emission 
cores, indicative of a small degree of chromospheric activ- 
ity. This is known to cause e xcess scatter in the velocities 
(e.g., ISaar fc Donahuel ll997D. and is often referred to as 
"jitter" . While we cannot rule out other systematics that 
may be causing our internal errors to be underestimated, 
we attribute most of the excess scatter to activity. For 
our final fit we have added 4.2 ms~^ in quadrature to 
our formal errors, which yields a reduced near unity. 
The parameters of the fit are not significantly changed, 
and are listed in Table [3l The minimum companion mass 
is Mp sini = 0.625 ± 0.015 x [(A4 + Mp)/Mq]^/^ Mjup, 
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Fig. 2. — (a) Radial-velocity measurements for HAT-P-3 along 
with our orbital fit, shown as a function of orbital phase. The 
center-of-mass velocity has been subtracted, (b) Bisector spans 
(BS) for our 10 Keck spectra (template + 9 iodine exposures), 
computed as described in the text. The mean value has been sub- 
tracted, and the vertical scale is the same as in the top panel. 

where represents the mass of the star. 

We investigated the possibility that the radial veloc- 
ities we measured are the result of distortions in the 
line profiles due to contami nation f rom an u nresolved 
eclipsing binary (jSantos et al. 2002; Torres et al. 2005), 
instead of being due to true Doppler motion in response 
to a planetary companion. We cross-correlated each 
Keck spectrum against a synthetic template matching 
the properties of the star, and averaged the correlation 
functions over all orders blueward of the region affected 
by the iodine lines. From this representation of the aver- 
age spectral line profile we computed the mean bisectors, 
and as a measure of the line asymmetry we calculated the 
"bisector spans" as the velocity difference between points 
selected near the t op and bottom of the mean bisectors 
(jTorres et al.ll2005l ). If the velocities were the result of 
a blend with an eclipsing binary, we would expect the 
line bisectors to vary in phase with the photometric pe- 
r iod with an amplitude sim ilar to that of the velocities 
(Queloz et al. 2001; Mandus hev et al]|2005f ). Instead, we 
detect no variation in excess of the measurement uncer- 
tainties (see Figure m bottom panel). We conclude that 
the velocity variations are real and that the star is or- 
bited by a Jovian planet. 

5. STAR AND PLANET PARAMETERS 

The mass and radius of the planetary companion scale 
with those of the parent star. We have estimated the stel- 
lar properties by comparison with stella r evolutio n mod- 
els from the Yonsei-Yale (Y^) series bv iYi et al.l (^001). 
However, because a direct estimate of the distance is un- 
available (HAT-P-3 was not observed during the Hip- 
parcos mission), we have used as a proxy for luminosity 
(L*) the quantity a/ Ri, from the light cur ve fit, which is 
closely related to the stellar density (see ISozzetti et al.l 
|2007( ). As described by those authors, a/R^, is typically 
a better constraint on L^, than the spectroscopic value of 
log (7i, which has a relatively subtle effect on the line pro- 
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files and whose determination is th erefore more suscepti - 
ble to systematic errors. Following ISozzetti et al.l (|2007l ) 
we have determined the range of stellar masses and radii 
that are consistent with the measured values of T^g, 
[Fe/H], and a/i?* within their uncertainties. Based on 
previous experience, for this application we have adopted 
more conservative errors for the spectroscopic quanti- 
ties of 80 K for Tcff and 0.08 dex for [Fe/H]. We obtain 
Ah = O-diQto^oll M© and = 0.824l°:g^^ R©. The 
predicted surface gravity, logt;^ = 4.577tQ 045, agrees 
well with the spectroscopic estimate (see Table [2]). The 
nominal age we infer for HAT-P-3 is only 0.4 Gyr, al- 
though the uncertainty is very large due to the unevolved 
nature of this K dwarf. While this nominal age would 
make HAT-P-3 the youngest known transiting extraso- 
lar planet host star, we note that there is no sign of 
the Li I A6708 line in the spectra of HAT-P-3 down to 
a factor of ~4 lower than the Li level in stars of simi- 
lar spectral type in the Hyades cluster, which have an 
age of ~0.7 Gyr. This would suggest an age for HAT- 
P-3 greater than 0.7 Gyr. Further evidence for an older 
age is provided by the measured chromosphcric activity 
index of the star based on the Ca II H and K lines in 
the Keck spectra: logW^^ is —4.59 ±0.06, i ndicating an 
age of ~1.3 Gyr based on the cahbration by lNoves et al.1 
(fl98l . 

Combined with the results from the light curve and 
radial velocity modeling, the above stellar parameters 
yield a planet mass of Mp = 0.599 ± 0.026 Mjup and a 
radius of Rp — 0.890 ± 0.046 Rjup. The surface gravity 
of the planet, whi ch is independent of th e stellar pa- 
rameters (see, e.g.. ISouthworth et al.l 120071 ). is log^p = 
3.310 ± 0.066 (cgs), and is consistent with the P/\oggp 
relation shown by those authors. These and other planet 
properties are listed in Table [3] 

The absolute visual magnitude we estimate for the star 
using the stellar evolution models {My = 5.86 ± 0.20) 
combined with the apparent visu al magnitude {V = 
11.561 ± 0.067: iDroege et al.l[200l yields a distance of 
140 ± 13 pc, ignoring extinction. With this distance es- 
timate we investigated the possibility that the faint 10" 
companion to HAT-P-3 mentioned in ^is physically as- 
sociated. We find that it is not: its Ks-ha.nd brightness 
as measured by 2MASS is approximately 2 magnitudes 
fainter than predicted for a star with its J — Kg color, 
using the best-fit isochrone to HAT-P-3. It is therefore 
a background star. 

6. DISCUSSION AND CONCLUDING REMARKS 

The properties of HAT-P-3b place it among the smaller 
transiting planets discovered so far, and are remark- 
ably similar to tho se of the recently announced X0-2b 
()Burke et al.ll2007D . not only in mass and radius (re- 
ported as 0.57 Mjup and 0.97 Rjup, respectively), but 
also in some of the characteristics of the parent star such 
as Teff and M^,. In the case of HAT-P-3b, with a mass 
of Mp = 0.599 Mjup, the radius measured here appears 
too small for a pure hydrogen-helium planet according to 
cu rrent theory. The mo dels for irradiated giant planets 
bv lFortnev et al.l (|2007f) indicate the radius could be ex- 
plained by the presence of a heavy-element core of some 
75 M© (~ i of the total Mp) for the nominal mass, or- 
bital semimajor axis, and age we infer, although uncer- 



tainties in these properties and in the planetary radius 
allow the range to be from a bout 50 to 100 M ^. HAT- 
P-3b thus joins HD 149026b (jSato et al.ll2005D as a hot 
Jupiter for which theory predicts a substantial fraction 
of heavy elements. Using the same models and the pub- 
lished properties of X0-2b, we find that X0-2b is also 
expected to have a 40 M© core representing about j 
of its total mass. These three he avy-element ric h plan- 
ets, along with others suggested bv lBurrows et al.l ([2007), 
provide important support for the core-accretion scenario 
of planet formation. HAT- P-3b appears consistent with 
the correlation proposed bv lGuillot et al.l ()2006f ) between 
the host star metallicity and the amount of heavy el- 
ements in the planet interior, although there seems to 
be at least on e exception to this relation (WASP-lb; 
IStempels et"all i2007'). Compared to HD 149026b, HAT- 
P-3b has a similar mass of heavy elements but consider- 
ably more hydrogen and helium, a feature that will need 
to be explained by planet formation models. 

As of this writing HAT-P-3b is the smallest extraso- 
lar transiting planet discove red photometrical ly. While 
HD 149026b and GJ 436b (jCillon et al.ll200l do have 
smaller radii, they were both originally discovered in 
Doppler survey s before transits were noticed. 

iGaudil (|2005l ) has argued that selection effects in signal- 
to-noise-limited transit surveys tend to bias the size 
distribution of close-in giant extrasolar planets toward 
larger radii, since those objects have deeper transits and 
are therefore easier to detect. The fact that the recent 
discoveries of HAT-P-3b and X0-2b are on the small 
side of the radius distribution would seem to contradict 
the argument. However, the transit depth in both of 
these cases is not particularly small (~1.4%), a result 
of the fact that the parent stars are K dwarfs rather 
than earlier type stars. Thus, there may still be a bias, 
although relatively small. We note also that the proper- 
ties of HAT-P-3b are entirely consistent with the Monte 
Ca rlo simulation s of th e population of transiting planets 
bv iFressin et al.l (|2007D (see, e.g., their Fig. 17), as are 
other recently announced transiting planets. This sup- 
ports the idea that these objects are quite representative 
of the general population of giant planets, and are not 
significantly different from the Doppler planets. 
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